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The Bochkov-Kuzovlev nonlinear fluctuation-dissipation theorem is used to derive 
Narayanaswamy’s phenomenological theory of physical aging, in which this highly nonlinear 
phenomenon is described by a linear material-time convolution integral. A characteristic property 
of the Narayanaswamy aging description is material-time translational invariance, which is here 
taken as the basic assumption of the derivation. It is shown that only one possible definition 
of the material time obeys this invariance, namely the square of the distance travelled from a 
configuration of the system far back in time. The paper concludes with suggestions for computer 
simulations that test for consequences of material-time translational invariance. One of these is 
the “unique-triangles property” according to which any three points on the system’s path form 
a triangle such that two side lengths determine the third; this is equivalent to the well-known 
triangular relation for time-autocorrelation functions of aging spin glasses [L. F. Cugliandolo and 
J. Kurchan, J. Phys. A: Math. Gen. 27, 5749 (1994)]. The unique-triangles property implies 
a simple geometric interpretation of out-of-equilibrium time-autocorrelation functions, which 
extends to aging a previously proposed framework for such functions in equilibrium (J. C. Dyre, 
cond-mat/9712222 . 


I. INTRODUCTION 

Physical aging is the small, gradual change of material properties observed for some systems, e.g., polymers and 
glasses just below the glass-transition temperature [l|. In contrast to the aging of most real-life materials resulting 
from chemical reactions, physical aging is exclusively due to adjustments of molecular positions. During the last 
50 years there has been slow, but steady progress in the description and understanding of aging in experiments as 
well as model systems. For example, studies of spin glass aging in the 1990s led to the concept of an “effective” 
temperature that quantifies violations of the fluctuation-dissipation theorem describing linear response around a state 
of thermal equilibrium ill- There is still no good theoretical justification, however, of Narayanaswamy’s remarkable 
phenomenological physical-aging theory from 1971 Q, which has been used routinely in industry for decades [l|. 

This paper is an attempt to understand the microscopic origin of Narayanaswamy’s theory in which physical aging, 
which is strongly non-linear in the temperature history, is described by a linear convolution integral if one replaces the 
time integration variable by the so-called material time. The material time quantifies how fast processes take place 
in an aging system, reflecting the existence of an “inner clock”. In the words of Narayanaswamy’s seminal paper [^: 
“The intrinsic relation between force and flow (or cause and effect) is assumed to be linear, i.e., it is assumed that 
the only cause of nonlinearity in glass transition phenomena is the changing viscosity of the glass as the structure (or 
fictive temperature) changes.” 

The Narayanaswamy theory, which works best for small temperature variations [I|, raises a number of questions: 

1. How and why can the highly non-linear physical-aging phenomenon be described by a linear convolution integral? 

2. What is the material time? Can it be defined in terms of the system’s microscopic variables? 

3. Are all aging quantities subject to the same material time? 

4. Which physical quantity controls the material time? 

The first three questions are addressed in this paper, which derives the Narayanaswamy theory from an exact, non¬ 
linear fluctuation-dissipation theorem, assuming material-time translational invariance (Sec. IIIII) . The microscopic 
expression for the material time is derived in Sec. |IVl Sections 0 and |Vl] discuss consequences and give a few 
suggestions for numerical tests of the ideas developed; Sec. IVIII concludes the paper. 


Electronic address: dyre@ruc.dk 





2 


We do not claim originality of any arguments and results of this paper, which is just an attempt to illuminate the 
Narayanaswamy 1971 aging theory by connecting it to the 1981 Bochkov-Kuzovlev non-linear fluctuation-dissipation 
theorem and the 1994 spin-glass aging theory of Cugliandolo and Kurchan [^. - Much has happened in physics and 
chemistry since 1971, but apparently the theory of aging itself matures only like a slowly aging system. 


II. BACKGROUND 

This section summarizes some necessary background. 


A. The material-time concept and Narayanaswamy’s phenomenological aging theory 

Physical aging is usually studied in glasses and glass-forming liquids. For these systems the relaxation time is 
dramatically temperature dependent (sometimes increasing up to a factor of ten for a temperature decrease of just 
one percent). Even the simplest case, that of a temperature juim from a state of equilibrium of a highly viscous 
liquid, is nonlinear unless the teimerature change is very small [l|. Another complication is that aging is virtually 
always nonexponential in time [I|, Q • 

Interest in physical aging has throughout the years been stimulated by the practical iiMortance of being able to 
predict aging quantitatively, and there is a large literature on physical aging experiments [l|, lfil-ll.8l| . An excellent and 
still not outdated introduction to the phenomenology of aging is given in Scherer’s book from 1986 [ij. Experiments 
have studied physical aging by probing, e.g. , density [1, [ij i enthalpy a a, Young’s modulus [l^, gas permeability 

a , high-frequency mechanical moduli flA [l8j| , dc conductivity Q , frequency-dependent dielectric constant |l9l - 
beta relaxation loss-peak frequency [2lj . X-ray photon-correlation-spectroscopy probed structure , non-linear 
dielectric susceptibility [^, etc. In some cases one monitors aging following a temperature jump from a state of 
equilibrium, but more complicated temperature histories have also been studied extensively, in particular that of 
continuous cooling of a liquid into its glass state. 

To be specihc we shall assume that the property X is monitored when temperature is changed in a controlled 
fashion, and that the experiment ends in equilibrium at the temperature Tq. The deviation of X from its equilibrium 
value at Tq is denoted by AX{t); thus AX{t) —> 0 as < —> oo. Narayanaswamy’s idea, which was inspired by rheological 
models, was to introduce a so-called material time, ^ Q. This quantity may be thought of as the time measured on 
a clock with a clock rate that itself ages sm, reflecting the fact that processes in an aging system proceed with a 
speed that changes as the structure ages. Eollowing a down jump in temperature, for instance, the clock rate decreases 
gradually; it eventually becomes constant, but only when the system is very close to equilibrium. The material time 
is intuitively analogous to the proper time of relativity theory, which is the time measured on a clock moving with 
the observer. 

If the clock rate is denoted by 7 (t), the material time f,{t) is defined by 


df = -f{t) dt. (1) 

In thermal equilibrium the clock rate is constant and the material time is proportional to the laboratory time, though 
with a proportionality constant K that is highly temperature dependent, i.e., ■feq(t) = K{T)t + Const. For an aging 
system, the clock rate is generally a complicated function of the system’s thermal history after it fell out of equilibrium. 

Narayanaswamy suggested that for the temperature variation T{t) = To + AT{t), if rewritten as a function of f the 
variation AX{t) has an “instantaneous” contribution, CAT{t) = CAT {fit)), and a relaxing contribution given by a 
linear material-time convolution integral [H, 0|i i-C-i 

dAT 

AA(6) = C'AT(6) - / Mx{f2 - fi) ^(6) df, . (2) 

J-oo dfi 

As an example, after a temperature jump at t = 0 from equilibrium at temperature Tq + AT to the temperature Tq 
one has AX{t) = ATMx{f{t)) in which f{f) is found by integrating Eq. ([T]). The crucial point is that the function 
Mx{f) is the same for all temperature jumps whereas the material time f{t) depends on AT (usually strongly). 

The Narayanaswamy formalism represented a great step forward compared to that time’s use of non-linear differ¬ 
ential equations for describing aging. Such equations cannot reproduce memory effects like the crossover (Kovacs) 
effect [l|, and they do not reflect the deep insight that aging can be described by a linear theory. Recently, the 
use of differential equations for describing aging has again been investigated, but now limited to temperature jumps. 
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Thus Kolvin and Bouchbinder have proposed a simple temperature-jump differential equation that fits data well [^ . 
and Hecksher et al. showed that for aging following temperature jumps the Narayanaswamy formalism is equivalent 
to a differential equation M- The latter differential equation involves both the normalized and the unnormalized 
relaxation function, and this is necessary to ensure consistency with Eq. ©• This equation results in a simple expres¬ 
sion for calculating one nonlinear relaxation function from another, without any fitting to analytical functions as is 
usually done in applications of the Narayanaswamy formalism. The resulting formalism was confirmed in experiments 
monitoring dielectric relaxation and the high-frequency shear modulus during aging [^ . 

In the Narayanaswamy theory the clock rate of the aging system’s “inner clock” is a global variable. An inter¬ 
esting extension of this formalism is due to Castillo and co-workers [26l - l^ who assume that aging is controlled 
by a space-dependent clock rate. This idea, which makes perfect sense in view of the existence of dynamical hetero¬ 
geneities in glass-forming liquids and glasses [1^ , has been verified numerically for different model systems [26l - [2^ . 

Nevertheless, below we do not follow this promising idea, but focus on the simpler case of a “global average” clock. 


B. Material-time translational invariance (MTTI) 

Equation ([2]) has the appearance of standard linear-response theory, though with the time t replaced by the material 
time f . A characteristic property of linear-response theory is time-translational invariance, which expresses the fact 
fundamental to all science that repeating an experiment at a later time leads to the same result. An implicit 
assumption of Narayanaswamy’s formalism is that translational invariance applies for the material time during aging 
as reflected in the fact the memory function in Eq. ([2]) is a function of the difference ^2 ~ Ci- This will be referred to 
as material-time translational invariance (MTTI). 

Translated into standard statistical-mechanical language, MTTI is the property that even during aging, time cor¬ 
relation functions depend only on the difference of the material times involved, i.e., that for any two quantities A and 
B a function Fab{0 exists such that 


{A{h)B{t2)) = FAB{at2)-ati)). (3) 

The angular brackets denote an ensemble average (for the aging system). In contrast to the equilibrium situation, 
the left-hand side is generally not a function of ^2 — • It is understood here and henceforth that A and B both have 

zero average - if that is not the case, the average value should be subtracted at each time t. 

Of particular interest is MTTI for the displacement vector that specifies the system’s path in configuration space, 
the 3A^-dimensional vector R = (ri, ...rjv) where is the position of the zth particle. One way to investigate a system’s 
dynamics is to evaluate the mean-square displacement from one time ti to a later time t 2 , ((R(t 2 ) ~ II(^i)) )■ This 
quantity, of course, refers to the total displacement, which impli es g oing beyond the boundaries of standard computer 
simulations performed with periodic boundary conditions With this clarification in mind, MTTI for the 

mean-square displacement means that it depends only on the difference in material times involved, i.e.. 


((R(t2) - R(ti))') = NGi^ih) - ati)) . (4) 

C. The triangular relation 

More than 20 years ago the study of aging in spin glasses resulted in a deep insight into the nature of aging in 
such discrete models Q. Interestingly, this work seems to have been carried out with little interaction with the 
experimental oxide glass and polymer aging communities, which had developed over long time 3,01. 

For spin glass aging there is time-scale separation, just as for the aging of ordinary glasses. On short time scales 
an analog of vibrational thermal equilibrium exists, whereas aging of a spin glass’ structural degrees of freedom takes 
place on much longer time scales. This led to the introduction of the concept of an effec tive temperature iim 
that is conceptually similar to Tool’s intuitive fictive temperature concept from 1946 (341) . but founded in rigorous 
statistical-mechanical theory. 

Moreover, the spin glass works resulted in the understanding that ’’the proper measure of time is the correlation 
itself, not the laboratory clock” [ 3 ^ . This was formalized into the triangular relation which, if C{t,t') is a time- 
autocorrelation function, i.e., an average for the aging system of the form C{t,t') = {A{t)A{t')), is the following 
equation valid on long time scales (ti <t 2 <tffj 0 
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C{h,t3) = f{C{h,h),C{t2,h)). (5) 

The idea that time is measured via the autocorrelation function itself UM implies time-reparametrization invariance, 
see, e.g.. Refs. [H, HI, H^, as well as the earlier Refs, [s^ andIsM 

As detailed shortly, the triangular relation is equivalent to MTTI, but first a few notes: 

1. Equation ([S]) applies in thermal equilibrium for any time-autocorrelation function that is a sum of decaying 
exponentials. In that case, the value of the time autocorrelation function C(ti, ^ 2 ) uniquely determines the time 
difference t 2 — ti, and C{t 2 ,t 3 ) likewise determines — ^ 2 , which means that knowledge of both C{ti,t 2 ) and 
C(t 2 ,t 3 ) implies that t^ — ti = (fa — 12 ) -f (t 2 — ti) and thus C(ti,fa) are known, which is the triangular relation 

Eq. (0. 

2. For the case of equilibrium, if time-temperature superposition (TTS) applies, the function / is the same when 
temperature is changed (for the normalized autocorrelation functions). Conversely, if Eq. ([5]) applies for aging 
over a certain range of temperatures, this implies TTS. 

3. The triangular relation applies whenever the system in question obeys so-called dynamic ultrametricity, i.e., 
when CiU t^) = mm{C{ti,t 2 ),C{t 2 ,t 3 )). This condition is believed to describe the aging of most or all spin 
glasses 

The triangular relation applies for any aging system that has an “inner clock” and a corresponding material time ^ 
obeying MTTI. In that case, if the function FcciO of Eq. dH) goes monotonically to zero as ^ goes to infinity, C{ti,t 2 ) 
determines the difference ^(^ 2 ) — Likewise, < 7 (^ 2 ,is) determines the difference — ^(^ 2 ). Thus C{ti,t 2 ) and 

Cit 2 ,t 3 ) determine ( 7 (^ 1 , < 3 ), which is the triangular relation. Note the analogy to the equilibrium argument. 

Conversely, if a system ages in such a way that its time-autocorrelation functions obey the triangular relation, one 
can define a material time ^(t) that satisfies MTTI. This nontrivial result was derived in Appendix B of Ref. We 
will not here repeat the derivation involving formal group theory, but briefly the argument is the following. Assuming 
again that C(ti,t 2 ) goes monotonically to zero as \t 2 — ti| —>■ 00 , one defines an algebraic operation * on the part 
of the positive real axis that corresponds to some value of the autocorrelation function: For x,y,z in this part of 
the real axis, we define x * y = f{x,y) where / is the function appearing in the triangular relation Eq. ([5]). It is 
straightforward to show from Eq. m that > 1 = is associative, i.e., that x * {y* z) = {x*y) * z. The commutative law also 
applies, i.e., x * y = y * x, which is not trivial - this signals an intriguing element of time reversibility during aging. 
From these facts one can show that functions h{t) and (j)(h) exist such that C{ti,t 2 ) = (t>{h{t 2 ) — [^. This is 

Eq. (O with h{t) = ^{t), i.e., the triangular relation leads to the existence of a material time and implies MTTI. 


III. DERIVING NARAYANASWAMY’S AGING THEORY FROM MTTI 

The fundamental puzzle of the Narayanaswamy formalism is how to arrive at a theory that describes aging by a 
linear material-time convolution integral. To show how this may come about we combine MTTI with Bochkov’s and 
Kuzovlev’s non-linear fluctuation-dissipation theorem from 1981 

It is convenient to adopt the language of the energy-bond formalism according to which any system interacts with 
its surroundings via one or more so-called energy bonds [dll - l^ . An energy bond is characterized by two variables, 
an effort and a flow variable, the product of which gives the energy per unit time transferred into the system from 
its surroundings. Effort could be, e.g., voltage drop or force, in which cases the corresponding flow variables are, 
respectively, electrical current and velocity. The formalism is general, and even a heat flow may be modelled - in 
that case energy is generalized into free energy, the effort is the temperature deviation from a reference temperature 
To, and the flow is the entropy current [46Ll47l|. The general nature of the energy-bond formalism means that for any 
physical quantity q one may define a flow / from / = dq/dt and let e be a fictitious field coupling linearly to q in the 
Hamiltonian. The thermal average effort and flow are assumed to be zero. 

We first consider the case of a single energy bond in which e(t) is the effort and f{t) the flow. If the effort is 
externally controlled, the standard linear fluctuation-dissipation (ED) theorem is the following expression for the 
average flow at time t 2 in the time-dependent external field e(t) 


{f{t 2 )) = /3o 


*2 


{f{t2)f{ti))oe{ti)dti. 


— 00 


( 6 ) 
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Here /3o = I/^bTq where Tq is the system’s temperature and (/(t 2 )/(H))o is the thermal equilibrium time- 
autocorrelation function, which by standard time-translational invariance is a function of the difference t 2 — ti- 
There is complete symmetry between e and / in the energy-bond formalism, and if the flow / is externally controlled 
the FD theorem is 


(e(t 2 )) =/3o / {e{t2)e{ti))Q f{ti)dti . (7) 

J —CO 

In 1981 Bochkov and Kuzovlev derived a general, non-linear response theory [d^. It is assumed that the system 
was in equilibrium at the temperature Tq far back in time and ends up in equilibrium at the same temperature in the 
distant future. If non-equilibrium cumulant averages are denoted by (...)c, the following exact relation (Eq. (2.18) in 
Ref. 113 in the energy-bond notation) applies: 

{f{t)) = X! ^ — ~f dti... j dtk{f{t)f{ti)...f{tk))ce{ti)....e{tk) ■ (8) 

To first order in e all but the first term on the right-hand side may be ignored and the two-point cumulant average 
may be replaced by the equilibrium product average {/(t 2 )/(^i))oi thus reducing Eq. ([S]) to the standard ED theorem 

Eq. dH). 

The derivation of Eq. ([5]) is based on Bochkov and Kuzovlev’s theorem [ 13 ] according to which the ratio of 
probabilities of a given path in time R(t) and the corresponding time-reversed path R(—<) is the ratio of the Boltzmann 
probabilities of the starting and ending points of the path (considered over some finite, large time interval). During 
the 1990s this and related results became known as the fluctuation theorem (see, e.g.. Refs. to the best of my 

knowledge its first general derivation was given in Refs. |13 and [ 13 , covering both the classical Newtonian, stochastic, 
and quantum-mechanical cases. 

Equation ([51) is not the simple higher-order generalization of the linear FD theorem in which the right-hand-side 
averages are equilibrium averages. In fact, as shown some time ago by Stratonovich [^ . there is no general nonlinear 
FD theorem determining the nonlinear response from equilibrium averages (40j| . The occurrence of non-equilibrium 
averages on the right-hand side makes Eq. ([8]) appear to be fairly useless. We proceed to show how one can, however, 
based on MTTI use Eq. ([5]) to arrive at the Naryanaswamy expression that is linear in the material time. 

The starting point is the already mentioned fact that Naryanaswamy’s theory works best for small temperature 
variations [l|. Eor larger jumps it breaks down; it is not known whether this is because multiple relaxation mechanisms 
must be taken into account, or whether there is a fundamental break down of the material-time idea (see, e.g., Ref. 
[53 and references therein). In any case, it is known from experiments that there is an aging regime described by Eq. 
(ED in which the temperature changes involved are fairly small (typically a few percent), but large enough that aging 
is strongly non-linear (e.g., with changes of a factor of ten or more in the characteristic relaxation time between up 
and down jumps to the same temperature 

Standard aging experiments are controlled by temperature, in which cases the effort e is the temperature difference 
to a reference temperature Tq. Eor the regime of relatively small temperature variations it makes good sense to keep 
only the lowest-order term on the right-hand side of Eq. ([5]). Unfortunately, as in most of physics, it is not possible 
to predict from Taylor expansions how good is the approximation of just keeping the first order term. This means 
that we cannot estimate how small the temperature jump must be for thisapproximation to be good. 

Ignoring the higher-order terms of Eq. ([5]) leads to 


ifih)) = Pq f {f{t2)f{ti))ce{ti)dti . (9) 

—00 

We have ignored the fact that the cumulant averages of Eqs. m and (ini) are themselves in general highly nonlinear 
functions of the effort variable e. The procedure of ignoring this is justified below a posteriori from MTTI and a 
change of integration variable to the material time. To consistently ignore higher-order terms one must also ignore 
the difference between the cumulant average (/(t 2 )/(U))c and the product average (/(t 2 )/(U)) (the difference is 
( 7 (^ 2 ))(/(U)) and of second order in e if (/(t))o = 0). This leads to 


*2 


{f{t 2 )f{ti))e{ti)dti. 


= /3o 


— 00 


( 10 ) 
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Equation (flOl) has the appearance of the standard linear FD theorem Eq. but since (/(t 2 )/(ii)) is evaluated 
along the actual system path, this quantity is as mentioned in general a highly nonlinear function of the effort history 
e{t) and totally different from any equilibrium time-autocorrelation function. The point of MTTI used below is that 
changing the integration variable to the material time eliminates this effort history dependence. 

We rewrite the flow variable / as a time derivative of a generalized “charge” q, 


f = 


dq 

dt 


and change the integration variable in Eq. (ITU)) from the time t to the material time ^ = ^{t). This leads to 


( 11 ) 





e(a)da 


( 12 ) 


By MTTI (Eq. ([3])) the autocorrelation function {q{^ 2 )q{^i)) is a function of ^2 — Cii i-^-, 


(9(a)9(6)) = i",,(6-ei)- (13) 

The crucial point is that all dependence on e(t) on the left-hand side comes from the material time’s effort dependence, 
whereas the function Fqq is independent of e{t). 

Equation (TTUl) implies 

and Eq. (HU thus becomes 

^(g(6)) = -/3o i""(6 -a)e(ei)rfa • (15) 

Replacing ^2 with ^ and integrating with respect to ^ from —00 to ^2 leads to (with Aq{^ 2 ) = 9 (^ 2 ) — 9eq(2o)) 


(Ag(6)) = -Po di f F" (^ - 6) e{S,^)di^ . (16) 

J —00 J —00 

Interchanging the order of integration we get since —00 < < ^ < ^2 

(Ag(6)) = -^0 dCi e(a) - i^)d^ . (17) 

J —00 J 

Utilizing that F^g(O) = 0, which follows from the fact that F(^) also describes the equilibrium linear response for 
which {d/dt){q{0)q{t))\t=o = 0 because of time-reversal invariance, we get 


(Ag(C2)) 



Cl) e(Ci) = Po 





e(Ci)dCi. 


When integrated partially, since Fqq{^) —J- 0 for ^ 00 , this leads to 


(18) 


(A(?(C2)) = Po 


(^Fqq{0)e{^2) 



Cl) ^(^1) 


(19) 


This is the sum of an “instantaneous” and a “relaxing” contribution. If the identifications q = X and e = AT are 
made, Eq. m is the Narayanaswamy expression Eq. ©■ This one cannot do a priori, however, because these 
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variables generally belong to different energy bonds, but the extension needed to derive Eq. ([2]) from Eq. dl) is 
straightforward, as we now proceed to show. 

Assuming that the effort e is the temperature variation AT - the usual aging situation - the corresponding flow 
variable is the entropy current, which is basicalfy the heat current into the system [d^. Equation ® generalizes 
straightaway to any number n of energy bonds [40|. The resulting equation involves a sum over all combinations of 
energy-bond indices, ji,...,jk, such that for i = each flow variable /ji(ti) in the right-hand-side cumulant of 

Eq. dHIl’s generalization is paired to the same-index effort variable ej^{ti). If there are n energy bonds, the leading 
term in the efforts is the following sum generalizing Eq. (nni) 


^ I't2 

(/(<2)) =/3oXI / (20) 

If one is interested in how the quantity X responds to a temperature variation, only one term in the above sum is of 
interest. This is the term coupling the X energy bond to the thermal energy bond whose generalized charge - the 
time-integrated entropy current - will be denoted by qs- In this way one arrives at an expression of the form Eq. 
(II9D with a different function T(^), but the same material time, thus affirmatively answering question 3 in Sec. HI 
The function F is given by F {^2 — ^i) = {qs{^i)X{^ 2 ))- Note that the above used identity T'(0) = 0 applies also in 
the general case because it follows from equilibrium time reversibility, {qs{ti)X{t 2 ))o = {X{ti)qs{t 2 ))o- 


IV. IDENTIFICATION OF THE MATERIAL TIME 

We proceed to show that there is only one possible material time ^{t) obeying the MTTI requirement. Consider the 
configuration-space path of an aging system, R(t). In the thermodynamic limit the relative fluctuations of distances 
go to zero. Because of this the ensemble-average symbols (...) on the left-hand side of Eq. may be removed, 
leading to the following formulation of MTTI for the displacement 


RI 2 = G{at2) - ati)), 


( 21 ) 


where we have introduced the notation 


,2 _ (R(t2)-R(ii))' 

■12 - AT 


( 22 ) 


The division by the number of particles N ensures a distance measure that is well defined in the thermodynamic limit 
{N —>■ 00 ). 

First, uniqueness of the material time is demonstrated. Equation (12111 applies, in particular, in equilibrium where 
7 (t) is a constant, implying ^(t) = Kt -|- Const. In thermal equilibrium, the generic behavior is that i ?22 is sublinear 
at short times, e.g., R 12 ^ (^2 — ti)^ with 0 < a: < 1, and linear in ^2 — ti at long times corresponding to ordinary 
diffusio n I55l - l57l| . The transition between the two regimes takes place around the system’s average (alpha) relaxation 
time iiHi; but we are not here interested in these details (s^ . The important thing is that the equilibrium mean- 
square displacement at long times is always proportional to time. From this and the fact that ^ = Kt Const, in 
equilibrium one concludes that 


G(e) oc C (e ^ 00 ). (23) 

This means that, if Ro is a reference configuration far back in time on the aging system’s trajectory, except for a 
multiplicative and an additive constant the only possible material time is defined by the squared distance from Rq to 
R(t) denoted by i?ot. 


m = Rlt, ( 24 ) 

If a dimensionless material time is wanted, one may multiply ^(t) with the particle number density to the power 2/3, 
but for simplicity we stick here to the above definition. 




Having established uniqueness of the material time, we proceed to show that the definition Eq. (l24l) is consistent. 
For this to be the case the difference in material time between two events on the aging system’s path in configuration 
space, 5 (^ 2 ) niust be independent of the choice of reference configuration. To show this we compare Rq 2 — Rqi 

with the analogous quantity using a different reference configuration, R(to)- Since R(ti) — R(tQ) = R(ti) — R(to) + AR 
in which AR = R(to) — R(io)j 


Rl,^ = + ARViV + 2(R(ti) - R(to)) • AR/N. (25) 

Likewise 

K 2 = Ro 2 + AR^N + 2{R{t2) - R{to)) ■ AR/N. (26) 

These equations lead to 

K 2 - Rl'i = RI 2 - Roi + 2(R(i2) - R(H)) • AR/N . (27) 

When to and tg are both far back in time, R(t 2 ) — R(ti) is uncorrelated with AR because in terms of the velocity 
V = R one has [R(t 2 ) — R(ii)] ■ AR = dsi J*° ds 2 V(si) • V(s 2 ), which goes to zero for to —>• —00 and to' ^ —00 
since velocities far apart in time are uncorrelated (even for an aging system). Thus the last term of Eq. (12711 vanishes 
in these limits, ensuring consistency of the material time definition Eq. (1241) . Note that the reference configuration 
need not, in fact, be one of the system’s distant past; any configuration far away may be selected as reference 
configuration. Note also that by differentiation of Eq. (24) one finds the following expression for the clock rate of Eq. 
(1): ^{t) = 2f_^V{t')-V{t)dt'/N. 


V. TWO CONSEQUENCES 

This section discusses two consequences of the formalism developed. 


A. The “unique-triangles property” 

The configuration-space path of a system in thermal equilibrium has an interesting geometric property. Consider 
three times, H < ^2 < to. Following Eq. (|22ll the corresponding distances between the configurations on the system’s 
path, R(ti), R{t2), and R(t 3 ), are denoted by R12, R13, and i? 23 - If one defines Req{t) as the distance that the 
equilibrium system travels over time t, i.e.. 


Req(t) = y^^^S, ( 28 ) 

we have R12 = Req(t2 — ti), R13 = Reqits — ti), and R23 = Req{t3 — t2)- Since Req{t) is an increasing function of t, 
this implies that the triangle formed by the points R(ti), R(t 2 ), and R(t 3 ) is unique in the following sense: If two 
of the triangle’s side lengths are known, the third one is also known. For instance, if R12 and R23 are known, i?i 3 is 
determined since i?i2 gives t2 — ti and R23 gives to — t2, from which to — ti = (to — t2) + (t2 — ti) and thus i?i3 may 
be deduced. A system for which any three points on its trajectory determine a unique triangle in the above sense will 
be referred to as obeying the “unique-triangles property”. 

It follows from Eq. (IMl) that the unique-triangles property also applies for an aging system. This is because the 
difference of two material times, ^(^ 2 ) — ^(^ 1)1 determines the distance between the corresponding configurations R(ti) 
andR(t 2 ), R 12 , and vice versa] likewise ^(t 3 )—^(t 2 ) determines i? 23 - Thus ^(ts)—^(ti) = (C(^ 3 )~C(i 2 )) + (C(^ 2 ) —C(^i)) 
on the one hand determines i?i3 and on the other hand is uniquely determined by R12 and i?23- ~ Note the close 
similarity to the triangular relation Eq. m i- 
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B. A geometric interpretation of time-autocorrelation functions 

It is a property of equilibrium dynamics that knowledge of the value of one time-autocorrelation function determines 
the time difference in question and thus all other time-autocorrelation functions. Via MTTI as expressed in Eqs. ® 
and o this property generalizes to aging systems; in particular it means that the value of {A{tx)A{t 2 )), which is a 
function of ^(^ 2 ) ~ is hr a one-to-one correspondence with the distance travelled, Ri 2 - Thus one can define a 

geometric autocorrelation function as the average of A{ti)A{t2) for all pairs of times of the aging system with the 
same distance i?i 2 , corresponding to the same difference in material time. We denote this geometric autocorrelation 
function by (A(O)A(i?)). MTTI leads to the identity 


{A{ti)A{t2)) 


{A{0)A{R)) 


R=Ri2 


(29) 


Since this applies also when the system stops aging and is very close to equilibrium, the geometric autocorrelation 
function (A(O)A(i?)) is identical to that characterizing equilibrium (Appendix). 

Equation (1291) generalizes a proposal for a geometric interpretation of equilibrium time-autocorrelation functions 
put forward some time ago [58l| . It was never published in a regular journal and is therefore briefly summarized in 
the Appendix. The idea is that the reason any time-autocorrelation function is small for large time separations is 
that the two relevant configurations are far from each other (we assume that A only depends on the system’s spatial 
coordinates, not the momenta). In glass-forming liquids all auto-correlation functions slow down in the same way as 
temperature is decreased, and moreover time-temperature superposition is often obeyed [H^. These facts are easily 
understood from the geometric interpretation of time-autocorrelation functions (Appendix), according to which the 
slowing down of Req(t) upon cooling controls all time-autocorrelation functions in the same way, compare Eq. (j29p . 
TTS arises if the geometric autocorrelation function is independent of temperature, which seems to be a reasonable 
assumption. 


VI. POSSIBLE NUMERICAL TESTS OF THE PROPOSED FRAMEWORK 

The MTTI assumption can be checked by computer simulations. To do this it is first necessary to identify a system 
that obeys the Narayanaswamy formalism. This may be challenging because computers are not yet able to simulate 
realistic aging situations, so temperature jumps larger than a few percent may be needed in simulations (^ . 

Once a suitable model system has been identified, the following tests can be performed: 

1. Predicting nonlinear aging from linear aging. The ultimate test of Eq. (I2|) is to investigate whether information 
from very small temperature-jump simulations, i.e., aging in the linear regime, is enough to predict aging 
following larger temperature jumps when the material time is defined by Eq. (1241) . 

2. The unique-triangles property for an aging system. This can be investigated in the same way the triangular 
relation is checked [^: First the system’s path in configuration space R(t) is traced out and stored. Different 
triplets of configurations, R(ti), R(t2), and R(t3), are then picked out and the corresponding distances i?i2, 
i?i 3 , and R 23 evaluated. The system obeys the unique-triangles property if for a given value of i?i 2 there is a 
one-to-one correspondence between R13 and R23. In practice, this may be checked by plotting R13 versus R23 
for a narrow range of i?i 2 values. 

3. The geometric ansatz for the time-autocorrelation functions. To check Eq. (1291) for an aging system one may 
proceed as follows. First evaluate the equilibrium geometric autocorrelation function (A(0)A(i?))o. This may be 
done by tracing out and storing the system’s path in configuration space R(f). For many pairs of configurations 
on this path, Ri and R 2 , the distance R 12 and product A(Ri)A(R 2 ) is evaluated; (A(0)A(R))o is then the 
average of these products with R = R 12 . After this, the same procedure is performed for an aging system to 
see whether the same geometric autocorrelation function applies. An equivalent test is to investigate whether 
AA 12 = A(R2) — A(Ri) is the same function of the distance R 12 in equilibrium and during aging, which can 
be done by plotting AA 12 versus i?i 2 and comparing the two plots. 

If simulations have also identified a model for which the Narayanaswamy theory does not apply, it would be 
interesting to see whether this system violates the unique-triangles property. This may illuminate whether that 
property is a good “thermometer” for which systems obey the Narayanaswamy theory and which do not. 
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VII. CONCLUDING REMARKS 

The purpose of this paper was to address the challenge of justifying from basic principles the phenomenological, 
but indisputably successful 44 years old Narayanaswamy description of aging. The above reasoning should not be 
regarded as an attempt to formulate a compelling theory for physical aging, but merely as suggesting one possible 
way to go about this subject. 

Taking material-time translational invariance (MTTI) as the starting point, we have seen that 

• MTTI is equivalent to the triangular relation discovered by Cugliandolo and Kurchan in 1994 from theoretical 
studies of spin glass aging Q. 

• MTTI in conjunction with the 1981 Bochkov-Kuzovlev exact nonlinear fluctuation-dissipation theorem leads to 
the Narayanaswamy aging theory for small temperature variations. 

• There is only one possible material time obeying MTTI, namely that defined by the distance squared to a 
configuration of the system’s distant past. 

• MTTI implies a geometric picture of time-autocorrelation functions for an aging system, according to which the 
time-autocorrelation is regarded as a spatial autocorrelation evaluated at the distance travelled in configuration 
space. 

An open question is the origin of the intriguing extension of Onsager reciprocity to aging systems’ time autocor¬ 
relation functions, i.e., the finding that in certain spin glass models {A{ti)B{t 2 )) = {B{ti)A{t 2 )) as discussed by 
Cugliandolo and Kurchan and by Franz and Virasoro [s^ . This applies trivially in equilibrium because of time- 
reversal invariance. Whenever Onsager reciprocity applies for an aging system, it indicates an element of reversibility 
during aging [H, HH, Given the analogous aging behavior of spin glasses and real glasses, it seems likely that 
Onsager reciprocity may also apply for the latter, meaning that aging for relatively small temperature variations has 
an element of reversibility. This is consistent with the fact that the algebraic operation * discussed in Sec. Ill Cl is 
commutative. 

In Sec. U we listed four questions relating to the Narayanaswamy formalism. The above developments addressed 
the first three of these. Thus we have shown I) how, by use of the Bochkov-Kuzovlev exact nonlinear fluctuation- 
dissipation theorem and MTTI, the highly nonlinear aging phenomenon may be reduced to a linear material-time 
convolution integral; 2) the material time was identified as the distance to a reference configuration far away; 3) 
whether all quantities age following the same material time is answered with a yes since the material-time definition 
of Eq. (l24ll is unique. 

The remaining question from Sec. Uis: what controls the material time? The traditional answer to this is the Active 
temperature [j, [^ . which by definition quantifies the structure in such a way that in Tj = T in equilibrium. 
The Active-temperature concept is hand waving, however, and it would be nice to have a microscopic understanding 
of what controls the material time’s clock rate jit) of Eq. (IT|) [3, [13, m. This amounts to solving one of the 
deepest problems in glass science: what controls the equilibrium relaxation time’s temperature dependence Aging 
studies may contribute to solving this problem by providing information beyond that obtainable from linear-response 
experiments. Eor instance, an aging experiment could in principle determine whether following a temperature jump 
there is an instantaneous change of the activation energy of the clock rate ■jlt) as predicted, e.g., by the shoving model 
iElsl, or whether there is no such instantaneous change, which is the expectation, e.g., from the Adam-Gibbs 
conAgurational-entropy model . Another example of aging studies with consequences for the general understanding 
of viscous liquid dynamics is the suggestion that the potential energy controls the material time and thus, in particular, 
glass-forming liquids’ equilibrium average relaxation time (^ . This interesting idea connects to and includes the 
material time of rheology [63, [b^, which was Narayanaswamy’s source of inspiration [^. 
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Appendix: Geometric theory of equilibrium time-autocorrelation fnnctions 

As in the main paper, the system’s path in conAguration space as a function of time is denoted by R(t) and we 
deAne the distance i?i 2 between two conAgurations R(ti) and R(t2) by Eq. (0^ . As mentioned, Auctuations become 
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insignificant in the N ^ oo limit, and consequently the distance between the system’s configurations at two times is a 
unique number. This simplification applies in equilibrium, as well as for an aging system. In thermal equilibrium the 
displacement of the system during time t is denoted by Req{t) (Eq. (12^ 1. Note that this quantity is experimentally 
accessible via the incoherent intermediate scattering function. 

Some time ago I proposed a geometric view of thermal-equilibrium time-autocorrelation functions The back¬ 
ground was the facts that time-temperature superposition (TTS) often applies for linear-response functions of glass¬ 
forming liquids and that, while the average relaxation time t{T) is always strongly temperature dependent, t(T) of 
different linear-response functions usually varies with temperature in exactly the same way [^ . A simple way to 
understand these facts is that there is a temperature-independent geometric autocorrelation function, 


(A(O)A(i?)). (30) 

by which is meant the quantity defined by averaging over all pairs of points along the system’s path in time separated 
by the distance R = Req{t), corresponding to the time interval t = t2 — ti between two configurations, R(ti) and 
R(t 2 )- This distance is unique for given t. Req{t) is an increasing function of t and thus in a one-to-one correspondence 
with t. In terms of the geometric autocorrelation function an equilibrium time-autocorrelation function is given 
by 


(A(0)A(t))o 


(A(O)A(R)) 


R=Rt,q{t) 


(31) 


As it stands, Ea. dSTl) is a tautology because the geometric average (A(O)A(i?)) is defined to make Eq. (I5T|) apply at 
the thermodynamic state point in question. 

One of the consequences of Ref. is that TTS finds a natural explanation if what happens when temperature is 
changed is simply the following: the geometric autocorrelation function {A(0)A(i?)) is unchanged, and Req{t) is also 
unchanged except for an overall scaling of time reflecting the slowing down upon cooling. This geometric “explanation” 
of TTS admittedly presupposes TTS for the mean-square displacement, but there is considerable evidence for this 
from theoretical and numerical studies of hopping in highly disordered landscapes . 

In Ref. [Hwe proposed an ansatz for calculating the geometric autocorrelation function (A(O)A(i?)) via a “double- 
canonical” statistical-mechanical average, and the theory was validated by simulations of a one-dimensional double¬ 
potential model. 

A simple example of Eq. m is the case of a Gaussian geometric autocorrelation function, (A(O)A(i?)) oc 
exp(—i?^/2i?Q), and linear diffusion, i?eq(0 = 6Dt, where D is the diffusion constant. This leads to an exponentially 
decaying time-autocorrelation function, the so-called Debye relaxation that is the simplest linear-response situation. 
A more realistic case also assumes a Gaussian geometric autocorrelation function, but combines this with the fact 
that the short-time mean-square displacement is usually subdiffusive . This leads to a high-frequency (short-time) 
non-Debye behavior. In experiments the relaxation is often Debye-like on the low-frequency side of the loss peak 
- at least for non-polymeric systems - which implies that a Gaussian decay of the geometric autocorrelation function 
must apply at long distances. 

The assumption that the equilibrium geometric autocorrelation function (A(O)A(i?)) is temperature independent is 
a quantitative expression of the physical idea that at different temperatures the system paths are very similar, except 
for the fact that the rate of motion has changed [l^. Whenever this applies, it is natural to expect that the same 
geometric autocorrelation function describes a system aging at these temperatures. 
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